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1.  INTRODUCTION 

A  continuing  challenge  in  the  design  of  new  energetic  materials  (e.g. 
explosives  and  propellants)  is  to  reconcile  the  necessary  metastability  with  the 
desired  insensitivity.  Simply  put,  they  should  be  sufficiently  unstable  to  be 
capable  of  suddenly  releasing  a  great  deal  of  energy,  but  stable  enough  that 
this  not  happen  before  it  is  wanted!  Since  these  two  objectives  are  intrinsically 
somewhat  contradictory,  the  realistic  goal  is  an  optimum  compromise,  which 
maximizes  energetic  performance  while  minimizing  sensitivity  to  unintended 
stimuli. 

Some  of  the  data  relevant  to  the  first  issue  can  now  be  predicted  with 
generally  satisfactory  accuracy;  this  includes  density  [1-3],  heat  of  formation 
[4-6],  detonation  pressure  and  velocity  [7,8]  and  specific  impulse  [9].  These 
properties  permit  meaningful  assessments  of  the  potential  level  of  energetic 
performance.  The  prediction  of  sensitivity,  however,  continues  to  be  an  area  of 
considerable  activity. 

A  key  point,  with  regard  to  both  issues,  is  the  decomposition  process  of  the 
compound:  What  are  its  energetics  and  how  readily  does  it  occur?  Our 
emphasis  in  this  chapter  shall  be  upon  factors  that  influence  the  ease  with 
which  decomposition  can  be  initiated  by  unwanted  external  stimuli,  i.e. 
sensitivity.  These  stimuli  may  be  of  various  t3T)es,  including  impact,  shock, 
friction,  heat  and  electrostatic  charge  [8].  Relative  vulnerabilities  to  these 
different  effects  need  not  be  the  same;  for  example,  the  onset  temperatures  for 
the  thermal  decomposition  of  TNT  (1)  and  HMX  (2)  are  quite  similar,  but  the 
latter  is  much  more  likely  to  undergo  detonation  upon  impact  [8].  It  has  been 
shown,  however,  that  there  is  a  general  correlation  between  impact  and  shock 
sensitivities  [10],  which  are  the  ones  upon  which  we  will  focus. 

Impact  and  shock  sensitivities  depend  upon  a  variety  of  factors:  chemical, 
structural  and  physical.  We  shall  limit  our  discussion  to  the  effects  of 
molecular  structure  (i.e.  chemical  composition  and  molecular  geometry). 


Thus  we  will  not  directly  address,  for  example,  the  mechanisms  of  “hot  spot” 
formation  in  the  crystal  [11-13],  or  the  role  of  particle  size  [10,12,14]. 


Impact  sensitivity  is  taken  to  be  inversely  proportional  to  the  height  (hso) 
from  which  a  given  weight  falling  upon  the  compound  has  a  50%  probability  of 
producing  an  explosion.  Shock  sensitivity  is  directly  proportional  to  the 
maximum  gap  width  through  which  a  standard  shock  wave  has  a  50% 
probability  of  causing  an  explosion.  It  is  important  to  recognize  that  both  types 
of  results  are  very  dependent  upon  the  specific  physical  conditions  of  the 
processes,  and  reproducibility  can  be  a  problem  [8,12,15,16]. 

2.  IMPACT/SHOCK  SENSITIVITY  AND  MOLECULAR  STRUCTURE: 

SOME  BACKGROUND 

2.1.  Structure-Sensitiviiy  Relationships 

The  effects  of  molecular  structure  upon  impact/shock  sensitivity  have  been 
analyzed  and  reviewed  in  some  detail  on  a  number  of  occasions  [8,10,12,17]. 
Over  the  years,  there  have  been  frequent  efforts  to  relate  the  experimental 
impact  or  shock  sensitivities  of  groups  of  compounds  of  a  given  type  (e.g. 
trinitroaromatics)  to  some  molecular  quantity  or  quantities.  Molecular 
stoichiometry  has  proven  to  be  remarkably  effective  in  this  context  [10,15,16,18- 
20],  most  notably  Kamlet’s  oxidant  balance  formula  [15,16,18],  which  is 
essentially  a  measure  of  the  oxygen  that  is  present  relative  to  what  is  needed  to 
convert  all  hydrogens  to  H2O  and  carbons  to  CO.  In  general,  the  larger  is  the 
oxidant  balance,  the  greater  is  the  impact  sensitivity. 

In  his  extensive  and  very  significant  analyses  of  the  impact  sensitivities  of 
energetic  compounds,  Kamlet  emphasized  the  roles  of  C-NO2  and  N-NO2 
bonds  as  “trigger  linkages”,  the  cleavage  of  which  is  often  a  key  step  in 
decomposition  processes  [15,18].  Considerable  experimental  evidence  has 
accumulated  in  support  of  this  view  [21-35].  Kamlet  also  argued,  with  regard 
to  C-NO2  compounds,  that  rotation  around  the  C-NO2  bond  has  a 
desensitizing  effect  since  it  reduces  the  amount  of  externally-provided  energy 


that  can  go  into  a  C-NO2  stretching  vibrational  mode  and  promote  bond¬ 
breaking  [18].  Thus  any  steric  hindrance  of  rotation  can  be  expected  to 
increase  sensitivity.  In  this  context,  it  is  relevant  to  note  recent  studies 
indicating  that  nitro  groups  have  an  enhanced  capacity  for  localizing 
transferred  vibrational  energy  [36,37]. 

Reflecting  this  emphasis  upon  the  C-NO2  and  N— NO2  bonds,  several 
studies  have  sought  to  correlate  impact  and  shock  sensitivities  with  properties 
of  these  bonds,  particularly  various  measures  of  their  stabilities.  These  have 
included  the  lengths  of  the  bonds  [38,39],  the  electrostatic  potentials  at  their 
midpoints  [38,40,41],  and  their  polarities  in  excited  states  of  the  molecules  [42]. 
Taking  a  different  approach,  Kohno  et  al  have  argued  that  the  impact 
sensitivities  of  nitramines  can  be  related  to  the  differences  between  the  N-NO2 
distances  in  the  gas  phase  and  in  the  crystal  [43,44]. 

2.2.  Some  Specific  Decomposition  Pathways 

The  various  structure-sensitivity  relationships  mentioned  above  typically 
treat  different  classes  of  compounds  separately.  Nitramines,  for  example,  are 
not  expected  to  fit  on  the  same  correlations  as  nitroaromatics  or 
nitroheterocycles.  Even  within  a  given  class,  however,  a  variety  of 
decomposition  mechanisms  may  be  operative.  Indeed,  Kamlet  and  Adolph 
found  it  necessary  to  establish  two  oxidant  balance  correlations  for 
nitroaromatics  [16],  one  being  only  for  molecules  having  a  CH-containing 
substituent  alpha  to  a  nitro  group,  3: 
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Evidently  the  a-H  is  likely  to  be  involved  in  the  early  stages  of  the 
decomposition  process  [16],  perhaps  moving  to  an  adjoining  nitro  group  to 
form  a  nitronic  acid  tautomer,  4  [35,45,46].  These  are  known  to  be  reactive  and 
unstable  [47].  The  transfer  or  loss  of  a  proton  to  yield  a  nitronic  acid  or  a 
nitronate  (aci)  anion  has  also  been  invoked  as  the  initial  step  in  the 
decompositions  of  other  energetic  molecules,  e.g.  picric  acid  [35,48]  and 
amine-sensitized  nitromethane  [49-51]. 

A  particularly  interesting  compound  is  TATB,  5,  which  shows  a 
remarkable  lack  of  sensitivity  [10,16].  The  decomposition  of  TATB  proceeds,  in 
its  early  stages,  through  furazan  and  furoxan  intermediates  [52,53],  e.g.  6,  the 


formation  of  which  is  believed  not  to  involve  any  significant  net  release  of 
energy  [10,54,55].  This  presumably  means  that  the  progression  to  detonation  is 
less  rapid  than  if  these  early  steps  were  exothermic.  This  has  been  proposed 
as  an  explanation  of  the  insensitivity  of  TATB  [10,55].  As  a  contrast,  TNT  (1), 
which  is  much  more  sensitive  than  TATB,  decomposes  through  the 
intermediate  7  (among  others)  [16,56,57],  which  is  formed  exothermically  [54]. 


5  6  7 


A  structural  feature  that  is  frequently  associated  with  instability  is  the 
presence  of  several  linked  nitrogens  [58].  Depending  upon  the  molecular 
environment,  this  can  provide  a  relatively  facile  decomposition  pathway 
through  the  loss  of  N2.  Storm  et  al  used  this  reasoning  to  interpret  the 
observed  sensitivities  of  some  picryl  triazoles  [59],  and  we  have  shown 
computationally  that  the  high  sensitivity  of  the  triazole  8  can  be  explained  in 
an  analogous  manner  [60].  On  the  other  hand,  some  derivatives  of  the 
tetraazapentalenes  9  and  10  have  shown  a  surprising  degree  of  stability 
[58,61,62];  we  have  speculated  that  this  is  related  to  the  relatively  positive 
character  of  the  two  triply-coordinated  nitrogens  [63] . 
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It  should  be  apparent  from  even  this  very  brief  discussion  that  impact  and 
shock  sensitivities  can  depend  upon  a  variety  of  chemical  and  structural 
factors.  Any  generalizations  should  be  made  cautiously,  and  are  likely  to  be 
subject  to  qualifications  and  limitations. 


3.  RELATIONSHIPS  BETWEEN  IMPACT  SENSITIVITIES  AND 
MOLECULAR  SURFACE  ELECTROSTATIC  POTENTIALS 


3.1.  Analysis  and  Characterization  of  Surface  Potentials 

We  have  shown  in  earlier  work  that  it  is  possible  to  quantitatively  relate  a 
variety  of  liquid,  solid  and  solution  phase  properties  to  the  electrostatic 
potential  patterns  on  the  siarfaces  of  the  individual  molecules  [64-66].  Among 
these  properties  are  pKa,  boiling  points  and  critical  constants,  enthalpies  of 
fusion,  vaporization  and  sublimation,  solubilities,  partition  coefficients, 
diffusion  constants  and  viscosities.  For  these  purposes,  we  take  the  molecular 
surface  to  be  the  0.001  au  contour  of  the  molecular  electronic  density  p(r), 
following  the  suggestion  of  Bader  et  al  [67]. 

The  electrostatic  potential  that  is  produced  at  any  point  r  in  the  space 
arotind  a  molecule  by  its  nuclei  and  electrons  is  given  by  eq  (1): 


V(r)= 


y  Za  _  f  p(r')dr' 
f|RA-r|  J  |r'-r| 


(1) 


Za  is  the  charge  on  nucleus  A,  located  at  Ra-  The  sign  and  magnitude  of  V(r) 
are  the  net  result  of  the  positive  and  negative  contributions  of  the  nuclei  and 
electrons,  respectively,  at  the  point  r.  For  V(r)  computed  on  the  molecular 
surface,  we  determine  the  local  maxima  and  minima  (i.e.  the  most  positive 
and  negative  values.  Vs, max  and  Vs, min),  the  average  deviation  IT  and  the  total 
variance  a^t-  The  latter  two  are  defined  by  eqs.  (2)  and  (3): 


n  =  l|;|v(ri)-vs| 

^  ft 


(2) 


(3) 


—  —  1  “ 

Vs  is  the  average  potential  over  the  entire  surface:  Vg  =  —  ^Vfrj).  V+fr^) 

^i=l 

and  V“(rj)  are  the  positive  and  negative  values  of  V(r)  on  the  surface,  and  Vg 

and  Vg  are  their  averages:  Vg  =  — 

^ i=l  j=l 


We  view  IT  as  an  indicator  of  the  internal  charge  separation,  or  local 
polarity,  that  is  present  even  in  molecules  with  zero  dipole  moments.  It  has 
been  shown  to  correlate  with  dielectric  constants  [68]  and  with  an 
experimentally-hased  measure  of  polarity  [65].  The  total  variance,  ofot,  which 
is  the  sum  of  the  positive  and  negative  variances,  reflects  the  spread,  or  range 
of  values,  of  the  surface  electrostatic  potential.  It  is  particularly  sensitive  to 
the  positive  and  negative  extremes,  because  of  the  terms  being  squared.  We 
have  found  a^t  to  be  effective  as  a  measure  of  a  molecule’s  tendency  for 
noncovalent  interactions;  for  example,  it  enters  into  our  expressions  for 
boiling  point  [69],  enthalpy  of  vaporization  [64],  solubilities  [70,71],  etc. 

3.2.  Unsaturated  C-Nitro  Derivatives:  Nitroaromatics  and  Nitroheterocycles 

The  electrostatic  potential  of  a  ground  state  atom  is  positive  everywhere 
[72];  the  nuclear  term  in  eq.  (1)  dominates  over  that  of  the  dispersed  electrons. 
When  atoms  combine  to  form  a  molecule,  some  negative  region  or  regions 
normally  develop;  these  are  often  due  to  lone  pairs  on  the  more  electronegative 
atoms,  e.g.  nitrogen,  oxygen,  the  halogens,  etc.,  but  they  may  also  reflect  other 
factors,  such  as  unsaturated  or  strained  C-C  bonds  [73,74].  The  introduction 
of  the  strongly  electron- withdrawing  nitro  group  into  an  unsaturated 
molecule  generally  has  the  effect  of  eliminating  the  regions  of  negative 
electrostatic  potential  due  to  the  it  electrons.  For  example,  whereas  benzene 
and  aniline  have  extensive  negative  potentials  above  the  ring,  nitrobenzene  is 
positive  ever3rwhere  except  near  the  oxygens  [75].  Analogous  statements  apply 
to  the  regions  above  the  triple  bonds  of  acetylene  and  nitroacetylene  [76]. 

Two  classes  of  unsaturated  compounds  that  are  of  interest  in  the  present 
context,  as  energetic  materials,  are  trinitroaromatics  and  nitroheterocycles, 
e.g.  imidazoles  and  triazoles.  The  electrostatic  potential  on  the  molecular 
surface  of  1,3,5-trinitrobenzene,  11,  which  is  the  parent  molecule  for  the 
trinitroaromatics  that  we  shall  consider,  has  local  maxima  above  the  ring  and 
neeu*  each  of  the  C-NO2  bonds  [77].  This  pattern  is  modified  somewhat  when 
other  substituents  are  introduced,  but  its  basic  features  remain  even  under 
the  influence  of  electron-donating  groups,  e.g.  in  2,4-diamino-l,3,5- 
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nitroheterocycles,  such  as  l-inethyl-l,4-dinitro-l, 3, 4-triazole,  13,  have  surface 
potential  maxima  near  the  C-NO2  bonds  [77].  The  only  negative  surface 
regions  shown  by  any  of  these  three  molecules,  11  - 13,  are  associated  with  the 
nitro  oxygens. 

The  surface  potential  maxima  near  the  C-NO2  bonds  are  a  particularly 
interesting  feature  of  this  linkage.  Buildups  of  positive  electrostatic  potential 
above  C-NO2  bond  regions  have  been  observed  in  a  variety  of  t3rpes  of 
molecules  [75-84],  and  it  has  been  demonstrated  that  they  can  serve  as 
channels  for  nucleophilic  attack  [78,84]. 

In  a  recent  study  of  the  computed  electrostatic  potentials  and  related 
quantities  on  the  molecular  surfaces  of  13  trinitroaromatic  derivatives  [77],  we 
noted  two  rough  trends;  both  the  local  polarity  IT  and  the  surface  potential 
maximum  above  the  ring.  Vs, max.  tend  to  increase  as  the  impact  sensitivity 
increases.  Accordingly  we  used  a  statistical  analysis  program  [85]  to 
investigate  whether  an  acceptable  quantitative  relationship  could  be 
developed.  We  found  that  the  sensitivity  h.50  could  indeed  be  represented  in 
terms  of  FI  and  Vs,max(ring);  our  best  expression  was  [77], 

h5o  =  a  [n2Vs,max(ring)]“l  +  +  y  (4) 

In  eq.  (4),  a  >  0,  P  >  0  and  y  >  0.  The  linear  correlation  coefficient  is  0.989  and 
the  standard  deviation  is  14  cm,  for  experimental  impact  sensitivities  that 
range  from  36  cm  to  320  cm.  In  view  of  the  uncertainties  associated  with  the 
measured  values,  we  were  pleased  to  obtain  as  good  a  correlation  as  this. 

Our  success  with  the  trinitroaromatics  prompted  us  to  try  the  same 
approach  with  five  nitroheterocycles  (2  nitroimidazoles  and  3  nitrotriazoles) 
[77].  These  molecules  do  not  have  surface  potential  maxima  above  the  rings, 
so  we  used  for  each  one  the  most  positive  of  the  maxima  near  its  C— NO2 
bonds,  Vs,max(C-N02).  These  tend  to  increase  as  the  impact  sensitivity 
increases,  just  as  was  observed  for  the  Vs,max(ring)  values  of  the 
trinitroaromatics.  While  IT  does  not  show  a  parallel  trend  for  these 
nitroheterocycles,  nevertheless  a  satisfactory  representation  of  hso  was 
obtained  involving  both  Vs,max(C— NO2)  and  11  [77]: 

h50  =  a  [n2Vs,max(C-N02)]-l  +  P  (5) 

In  eq.  (5),  a  >  0  and  P  >  0.  The  linear  correlation  coefficient  is  0.986  and  the 
standard  deviation  is  19  cm,  for  experimental  sensitivities  between  35  cm  and 
291  cm.  There  is  a  notable  similarity  between  the  first  terms  on  the  right  sides 
of  eqs.  (4)  and  (5). 


A  possible  explanation  for  the  role  of  Vs,max(C-N02)  in  eq.  (5)  is  provided 
by  our  finding  that  it  shows  a  fair  correlation  with  the  computed  dissociation 
energy  of  the  C-NO2  bond,  for  a  total  of  eight  such  bonds  in  a  group  of  five 
nitroheterocycles  [86].  When  we  sought  to  generalize  this  observation  to  a 
series  of  eleven  nitroalkanes,  we  did  find  a  relationship  between  the  calculated 
C-NO2  dissociation  energy  and  Vs,max(C-N02),  which  now  also  includes  the 
molecular  surface  area  [87].  (The  fact  that  the  latter  was  not  needed  for  the 
nitroheterocycles  may  mean  simply  that  those  five  molecules  have  similar 
areas.)  Thus  it  may  be  that  Vs,max(C-N02)  in  eq.  (5)  is  an  indirect  measure  of 
the  strengths  of  the  C-NO2  bonds  in  the  nitroheterocycles;  this  would  be 
consistent  with  the  widespread  emphasis  upon  C-NO2  and  N-NO2  bonds  as 
trigger  linkages,  mentioned  earlier.  For  the  nitroaromatics,  it  may  be  that  the 
potential  maximum  above  the  ring  reflects  all  three  C-NO2  bonds,  due  to  their 
symmetrical  distribution. 

Understanding  the  function  of  n  in  eqs.  (4)  and  (5)  poses  a  greater 
challenge.  Eq.  (5)  and  the  very  rough  correlation  between  n  and  sensitivity 
among  the  nitroaromatics  seem  to  suggest  that  the  latter  increases  with 
internal  charge  separation,  which  presumably  counteracts  the  stabilizing 
effect  of  electronic  delocalization  in  these  unsaturated  molecules.  We  have 
advanced  this  argument  in  the  past  [77,87].  For  the  nitroheterocycles, 
however,  even  though  IT  does  appear  in  eq.  (5),  its  values  do  not  correlate  in 
any  obvious  manner  with  sensitivity  [77].  The  form  of  eq.  (4)  also  raises  a 
question:  The  contributions  of  the  two  variable  terms  on  the  right  side  are  both 
positive  and  of  the  same  order  of  magnitude  [77];  yet  11  appears  in  the 
numerator  of  one  and  the  denominator  of  the  other.  This  ambiguity 
concerning  the  role  of  11  led  us  to  reexamine  the  relationship  of  impact 
sensitivity  to  the  molecular  surface  electrostatic  potential. 

3.3.  Impact  Sensitivity  and  Surface  Potential  Imbalance 

The  negative  portions  of  a  molecular  surface  are  generally  the  smaller  part 
of  the  total  area,  but  they  are  frequently  relatively  strong;  thus  the  average 
negative  surface  potential,  Vs“,  is  typically  larger  in  magnitude  than  its 
positive  counterpart,  Vs^,  and  the  negative  variance  o?  is  greater  than  the 
positive,  a+.  This  can  be  seen,  for  a  representative  group  of  molecules,  in 
Table  1.  The  exceptions  tend  to  be  those  having  several  strongly  electron- 
attracting  constituents,  such  as  CF4  and  p -dinitrobenzene.  It  is  instructive  to 
compare  the  last  of  these  to  nitrobenzene;  the  second  nitro  group  actually 
decreases  the  magnitudes  of  both  Vs~  and  <t?,  because  the  polarizable 
electronic  charge  must  now  be  shared  between  two  electron-withdrawing 
substituents.  We  have  encountered  analogous  situations  on  a  number  of 
previous  occasions  [76,81,88].  One  consequence  is  that  the  value  of  n  tends  to 
level  off  as  the  number  of  strongly  electron-attracting  constituents  increases 


Table  1. 


Some  computed  molecular  surface  quantities.^ 


Molecule 

n 

Vs^ 

Vs 

a? 

cyclohexane 

2.2 

2.7 

-1.6 

2.5 

0.7 

C6H5CH3 

4.6 

4.2 

-5.2 

6.8 

11.1 

benzene 

4.9 

4.8 

-5.0 

7.1 

9.2 

(H5C2)20 

6.7 

5.9 

-9.4 

8.0 

129.8 

(HaOsCOH 

7.7 

6.2 

-12.3 

31.1 

182.7 

CF4 

8.3 

11.5 

-4.6 

66.9 

2.9 

pyridine 

8.5 

6.8 

-13.0 

18.5 

212.3 

CeHgOH 

8.6 

8.7 

-8.5 

63.8 

73.7 

CH3COCH3 

9.4 

7.7 

-18.8 

15.9 

159.8 

CH3COOCH3 

10.0 

7.4 

-16.7 

9.7 

129.2 

C2H5OH 

10.1 

8.1 

-13.7 

45.1 

182.4 

(H3C)2NCH0 

11.1 

9.3 

-17.0 

18.6 

158.8 

pyrimidine 

12.0 

9.9 

-17.4 

24.4 

163.0 

C6H5NO2 

12.3 

10.4 

-22.1 

16.7 

105.2 

CH3CONH2 

14.7 

12.5 

-20.0 

67.9 

139.9 

P-C6H4(N02)2 

16.5 

17.9 

-17.2 

29.8 

62.5 

H3CCN 

17.1 

15.5 

-22.2 

23.6 

167.8 

H2NCHO 

18.1 

17.1 

-20.0 

85.5 

233.6 

CH3NO2 

19.9 

19.4 

-21.4 

34.4 

81.7 

H2O 

21.6 

19.8 

-24.3 

85.7 

161.8 

^Data  are  taken  from  references  65  and  68.  n,  Vs^  and  Vs  are  in  kcal/mole; 
o+  and  al  are  in  (kcal/mole)^. 


[89];  thus,  in  going  sequentially  from  benzene  to  tetranitrobenzene,  the 
magnitudes  of  11  are  4.9,  12.3,  16.5,  19.5  and  21.4  kcal/mole.  It  is  actually 

larger  for  H2O,  21.6  kcal/mole,  than  for  tetranitrobenzene! 

Energetic  molecules  generally  have  several  strong  electron  attractors, 
such  as  nitro  groups  and  aza  nitrogens,  that  are  competing  for  the  polarizable 

charge.  Accordingly  they  tend  to  have  Vs  >  |Vs  |  and  a+  >  a_,  as  can  be 
seen  for  the  examples  in  Table  2.  This  is  in  marked  contrast  to  the  more 
typical  situation,  illustrated  in  Table  1.  For  the  former,  n  levels  off  in  the 
range  23  -  26  kcal/mole,  and  therefore  cannot  reflect  the  increasing  sensitivity 


Table  2. 

Measured  impact  sensitivities  and  computed  molecular  surface  quantities  for 
some  energetic  molecules.^ _ 

Molecule  1150  11  Vs^  Vs" 

nitroaromatics: 


X,Y.Z 


CHO,  H,  H 

36. 

19.1 

23.8 

-13.9 

152.3 

49.0 

OH,  H,  H 

87. 

20.9 

24.9 

-17.0 

174.3 

87.5 

H,  H,H 

100. 

19.5 

23.9 

-15.3 

109.0 

55.3 

CH3,  H,  H 

160. 

18.3 

22.1 

-15.2 

95.5 

51.7 

NH2,  H,  H 

177. 

18.6 

22.0 

-16.1 

97.9 

65.5 

NH2,  NH2,  H 

320. 

17.1 

19.4 

-16.8 

66.2 

75.2 

nitroheterocvcles : 

H 

1 

O2N.C.N.C.NO2 

\\  // 

C-N 

02N^ 

68. 

20.7 

27.4 

-12.9 

299.4 

43.7 

H 

1 

^N02 

w  // 

C-N 

O2N" 

105. 

23.3 

28.6 

-17.9 

312.1 

70.6 

H 

1 

NOo 

291. 

20.1 

23.4 

-16.4 

227.2 

79.7 

C'  C'  ^ 
\  n 

N-N 

/ 

H 

nitramines: 

NH-NO2 

H2C' 

NH-NO2 


13. 


23.4  29.6  -16.8  319.4  60.6 


NOo 
I  ^ 


HoC  CHo 

2|  j  2 

O^N'  g 
^2 

26. 

23.2 

27.8 

-18.5 

181.1 

57.2 

n  NO2 

C-N-CHg 

79. 

16.4 

18.3 

-14.4 

67.3 

47.7 

^C-N-CHo 

0  ' 

^  NO2 

CH3 

02N-NH-(CH2)2-N-N02 

114. 

22.6 

24.2 

-21.4 

164.0 

85.6 

N020 

1  ^  II 

H3C-N— CH 

320. 

17.4 

18.8 

-16.3 

51.2 

57.6 

^Impact  sensitivities  (in  cm)  are  from  reference  10;  the  computed  data  are 
from  reference  89.  fl,  Vs  and  Vg  are  in  kcal/mole;  a+  and  are  in 

(kcal/mole)2, 

that  accompanies  a  continuing  introduction  of  electron-withdrawing 
constituents. 

Our  analysis  led  eventually  to  the  conclusion  that  impact  sensitivity  can  be 
related  to  this  characteristic  imbalance  between  the  larger  Vg^  and  the 
smaller  |  Vg”  |  [89] .  There  are  of  course  a  variety  of  possible  measures  of  such 
imbalance,  including  the  quantity  v,  defined  by  eq.  (6),  which  we  have  used 


effectively  for  this  purpose  in  the  past  [64-66].  By  definition,  v  attains  a 
maximum  value  of  0.250  when  <t+  =  o?,  but  approaches  zero  when  »  a? 
or  a?  »  CJ+.  Among  other  possibilities  are  ratios  and  differences  of  Vg^  and 
I  Vg~  I  or  CT+  and  a?.  We  investigated  a  number  of  options  [89],  and  eventually 
settled  upon  eqs.  (7)  -  (9)  to  represent  the  impact  sensitivities  of  the  original  13 
nitroaromatics  and  5  nitroheterocycles,  as  well  as  a  group  of  8  nitramines. 
The  good  correlation  coefficients  and  relatively  low  standard  deviations  are 
gratifying,  given  the  uncertainties  in  the  data  bases. 


nitroaromatics: 


^50  -oc 


^  1 


l^^S.max  J 


+  P 


m  -m 

a,  P  >  0;  Y  <  0.  Correlation  coefficient  =  0.990;  standard  deviation  =  14  cm. 


+  Y 


(7) 


nitroheterocycles: 


^50  -ot 


Ve 


+ p 


(4 


+  Y 


(8) 


a,  P  >  0;  Y  <  0.  Correlation  coefficient  =  0.998;  standard  deviation  =  8  cm. 


nitramines: 


^50 


■  a 


^2  ^2 
cy+  <y_ 


«.2  „2 

—  c_ 


+  p 


<7+ 


+  Y 


(9) 


a  >  0;  P,Y  <  0.  Correlation  coefficient  =  0.997;  standard  deviation  =  9  cm. 


There  is  a  unifying  element  in  these  equations,  in  that  they  all  focus  upon 
the  relative  magnitudes  of  Vs^  and  |Vs~|  or  and  <y^.  The  additional 
dependence  upon  Vs, max  in  eq.  (7)  may  be  due  to  the  S3anmetrical  trinitro¬ 
benzene  framework  that  is  common  to  all  13  nitroaromatics.  In  eqs.  (8)  and 
(9),  the  contribution  of  the  first  variable  term  is  larger  than  that  of  the  second 
by  an  order  of  magnitude;  the  latter  can  be  viewed  as  a  correction  term. 

A  general  conclusion  that  seems  to  be  indicated  by  these  results  is  that  the 
anomalous  imbalance  between  the  strengths  of  the  positive  and  negative 
surface  potentials  is  at  least  S3nnptomatic  of  the  degree  of  instability  within 
these  three  classes  of  molecules.  It  should  be  noted  that  the  very  strong 
positive  potentials  produced  by  the  electron-withdrawing  components  do  not 
lead  to  as  high  levels  of  internal  charge  separation  (and  values  of  11)  as  might 
be  anticipated,  because  the  negative  potentials  are  uncharacteristically  weak. 


4.  SUMMARY 

We  have  presented  an  overview  of  various  attempts  to  relate  the  impact  and 
shock  sensitivities  of  energetic  materials  to  their  molecular  structures.  The 
objectives  of  such  efforts  are  to  better  understand  the  chemical  and  structural 
determinants  of  these  sensitivities,  and  to  develop  a  predictive  capability  to 
facilitate  the  evaluation  of  new  and  proposed  energetic  compoimds. 


Our  particular  emphasis  in  this  discussion  has  been  upon  the  relationship 
of  impact  sensitivities  to  the  electrostatic  potential  patterns  on  the  molecular 
surfaces.  The  current  status  of  our  analyses  is  represented  by  eqs.  (7)  -  (9). 
While  the  success  of  these  expressions  is  pleasing,  we  certainly  do  not  claim 
that  they  are  in  final  form.  They  reflect  small  data  bases  and  measurements 
with  a  relatively  high  level  of  uncertainty.  As  more  compounds  are  included, 
it  may  well  be  that  the  specific  formulations  given  in  eqs.  (7)  -  (9)  will  be 
modified;  i.e.  other  functions  of  Vs“  and  Vs"^  or  a?  and  aj  may  turn  out  to  be 
more  effective.  What  is  important  at  present,  however,  is  the  unifying  concept 
that  we  have  found  to  apply  to  all  three  of  these  classes  of  compounds,  namely 
that  their  impact  sensitivities  can  be  related  to  the  degree  of  imbalance 
between  their  typically  stronger  positive  surface  electrostatic  potentials  and 
weaker  negative  ones. 
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